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Abstract 



A one-parameter family of random variables, called the Discrete Action, is 
defined for a 2-dimensional Lorentzian spacetime of finite volume. The sin- 
gle parameter is a discreteness scale. The expectation value of this Discrete 
Action is calculated for various regions of 2D Minkowski spacetime, M^. 
When a causally convex region of is divided into subregions using null 
hues the mean of the Discrete Action is equal to the alternating sum of the 
numbers of vertices, edges and faces of the null tiling, up to corrections that 
tend to zero as the discreteness scale is taken to zero. This result is used to 
predict that the mean of the Discrete Action of the flat Lorentzian cylinder 
is zero up to corrections, which is verified. The "topological" character 
of the Discrete Action breaks down for causally convex regions of the flat 
trousers spacetime that contain the singularity and for non-causally convex 
rectangles. 
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1 The 2-dimensional action of a causal set 

The twin hypotheses that spacetime is fundamentally discrete and that, of all the 
structures of classical General Relativity, it is the causal structure of spacetime that 
will persist in the deep quantum regime gives rise to the idea that spacetime is a discrete 
order [H 121 [3] . Indeed the basic proposal of the causal set approach to quantum gravity 
is that the sum-over-histories for quantum gravity is a sum over discrete orders or causal 
sets. To define such a sum-over-histories, it will be necessary to give the amplitude for 
each causal set (or pair of causal sets if the path integration is conducted in Schwinger- 
Kel'dysh manner) and progress has recently been made on the question of what these 
amplitudes might be: a 2-dimensional action and a 4-dimensional action for a causal 
set have been proposed [1] and actions in 3, 5 and higher dimensions can also be defined 



Recall that a causal set {causet for short) is a locally finite partial order, i.e. it is a 
pair (C, :<) where C is a set and ^ is a relation on C which is reflexive {x ^ x), acyclic 
{y^x^y^y = x) and transitive (z^y^x^z^x). Local finiteness is the 
condition that the cardinality of any order interval is finite, where the (inclusive) order 
interval between a pair of elements y ^ x is defined to be I{x, y) '.= {z E C\y ^ z ^ x}. 
We call X the top element and y the bottom element of I{x,y). We write y -< x when 
y ^ X and y ^ x. We define n{x,y) := \I{x,y)\ and call a relation y -< x a. link if 
n{x,y) = 2. A chain is a totally ordered subset of C. 
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Sprinkling is a random process that produces a causet which is a discretisation of a 
d-dimensional, causal, Lorentzian manifold {Ai,g). It is a Poisson process of selecting 
points in {Ai,g), independently at random, with density p so that the expected number 
of points sprinkled in a region of spacetime volume V is pV. In quantum gravity we 
expect that the density is Planckian so that p = where / is of order the Planck 
length, but in this paper we treat p as a parameter to be varied. This process generates 
a causet whose elements are (identified with) the sprinkled points and whose order is 
that induced by the manifold's causal order restricted to the sprinkled points. If {M.,g) 
is of finite volume, the causet generated is almost surely finite and so the process defines 
a probability distribution Px,g,p on the set of finite causets {aka the set of finite partial 
orders). Henceforth, for ease of notation, we will drop the explicit reference to the 
metric g and refer, for example, to a spacetime as Ai and the probability distribution 
above as fM,p- 

Sprinkling is not a physical process. It plays a purely kinematical role and ex- 
presses the discrete-continuum correspondence: a causet C is well approximated by a 
Lorentzian manifold M. if it could have been generated, with relatively high probability, 
by sprinkling into Ai. In other words, C is well approximated by a manifold Ai if there 
exists an embedding i : C ^ M. such that (i) x, ?/ G C, ?/ ^ x iff iiy) G J^{i{x)) and 
(ii) the number of elements embedded in any sufficiently nice, large region of volume 
V is approximately pV . Strictly, this is a conjecture, the "Hauptvermutung" of the 
causal set approach, but it is supported by much evidence including the result that 
a distinguishing Lorentzian geometry is fully determined by its causal structure and 
spacetime volume measure O [71 [8] . 

We define the 2D action, S", of a finite causal set C to be [1] 

5[C] = iV - 2iVi + 4A/'2 - 2A/'3 (1.1) 

where N is the cardinality of C, and is the number of inclusive order intervals in C 
of cardinality m + 1. A^^i therefore is the number of links in C, N2 is the number of order 
intervals that are 3-chains (3 element chains) and A^4 is the number of order intervals 
that are 4-chains plus the number that are "diamonds" (with two mutually unrelated 
elements between the top and bottom elements). Note that is not the number 
of subcausets that are 3-chains but the number of order intervals that are 3-chains. 
The form of 5* as an alternating sum of (weighted) numbers of things is intriguingly 
reminiscent of certain topological indices. 



The action (1.1 ) defines an integer valued random variable, the Discrete Action Sm 



pi 
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for each finite volume spacetime Ai and density p via the sprinkhng process: Sm,p takes 
the value S[C] with probabihty P_A4,p(C). We also define the random variable Sm,n 
which takes the value S[C] with the probability that causet C arises in the process of 
selecting exactly elements uniformly at random - according to the spacetime volume 
measure - from Ai. We then have 



N 

= E '-^e-'^''{^M,N) (1.2) 



ipv) 

N=0 



where (■) denotes the expected value, V is the spacetime volume of Ai, and ^^]^| e~^^ 
is the probability that elements are selected in the Poisson process of sprinkling into 
At at density p. 

The Poisson distribution gives for the mean, 



{Sm,p) =pV- 2p2 [ d^y^/-g{y) [ d^x^/-g{x) 

Jm JMnJ+(y) 

(^1 - 2pV,y + ^{pV^yf^ e-^^- (1.3) 

where Vxy is the volume of the spacetime causal interval, [x,y] := J^{y) H J^(x), be- 
tween X and y and d is the dimension of A^. This can be understood thus: p d'^x^J —g{x) 
is the probability that an element is sprinkled in an elemental volume at x and similarly 
for y; e"''^="«, pK^e"''^^^ or \{pVxyY e~ P^^'y is the probability that there is no element, 
one element or two elements, respectively, sprinkled in [x,y\. 
Note that the double integration may be done in either order: 



{Sm,p) =pV- 2p2 / d^x^/-gix) [ d^y^-g^y) 
Jm JMnj-(x) 

(l-2pV;., + ^(pK.,)')e-^^-. (1.4) 



Indeed, the causet action (1.1) is invariant under reversal of the order relation on C, 
and so the Discrete Action (DA) for any spacetime {Ai,g) is equal to the DA of its 
t ime-or lent at ion-reverse . 

Sm,p is defined for any finite volume (causal) spacetime of any dimension so we can 
ask in what sense it is 2-dimensional. Each realisation of Sm,p is the action S[C] of 
some finite causet C and 

5[C] = 5^L(e.) (1.5) 

ei£C 
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where 



L{ei) = 1 - 2ni(ei) + 4n2(ei) - 2ns{ei) 



(1.6) 



and nm(ej) is the number of inclusive order intervals in C with cardinality m + 1 and 
with top element Cj. L(-) itself defines a random variable, Lj^^^py, for each spacetime 

, each point y E Ai and each p in the following way. Fix ?/ G , sprinkle into at 
density p and add an element at y to the sprinkled causet to form causet C which has 
a marked element, call it Cy. The value of 'LM,p,y is then L{ey) evaluated in C . If M. is 
2-dimensional, the mean of ^iM,p,y tends to j^R{y), where R{y) is the Ricci scalar, as p 
tends to infinity jl]. It approaches its limit when the discreteness length scale / := p~^ 
is much smaller than the curvature scale R~^. If Ai is not 2-dimensional, there is no 
apparent reason for Lj^^p^y to have anything to do with the continuum geometry M. 

Since L is thus related to the Ricci scalar when the causal set is a 2D sprinkling and 
S* is a sum of L(-) over the causal set, this implies that when Ai is 2-dimensional and 



terms arising from boundary effects. We will investigate this and in particular the na- 
ture of the boundary terms. In doing so we will be exploring whether the 2D Discrete 
Action is topological in character. The standard gravitational action for 2D Euclidean 
gravity, with its Einstein-Hilbert term and the (2D analogue of the) Gibbons-Hawking 
boundary term, is known to be a topological invariant, due to the Gauss-Bonnet the- 
orem. The Gauss-Bonnet Theorem has been extended to Lorentzian manifolds [91 [TO], 
so for ordinary (Lorentzian) 2D gravity, the action with an appropriate boundary term 
is also topological and a question arises: to what extent is the 2D causal set action 
topological? 



Consider a causal interval in 2D Minkowski spacetime, X := [p,q] C M^. For definite- 
ness consider the interval to have fixed volume (area), V. 

Following a conjecture of R. Sorkin, G. Brightwell proved that the mean {Sx,n) = 1, 
for any N ^ ^llj. This implies that the mean of Sjp is 



as p -> oo, {Sm,p) will tend to something that contains a term 




2 Intervals in 




= l-e-P^ (2.1) 

where ^^^j— e"''^ is the probability, in the Poisson process, that N elements are sprin- 
kled into X. 
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We use (1.3) to prove this result in a different way: 

(Sx) =pV- 2p2 f £y [ d^x p{pV,y) (2.2) 
Jx JxnJ+{y) 

where p{^) = (1 — 2^ + ^^^) exp(— ^) and we have suppressed the subscript p on the 
random variable Si^p. 

We use coordinates in which p and q lie on the time axis and q is at the origin. We 
consider null coordinates = -^{x'^ ~ x^), '^x = :^(^° + similarly for Uy, Vy. 

Then the interval is defined by m, f G [0, a] for a = y/V. 

pa pa pux pvx 

(Sx) = pV — 2 I dux / dvx I duy I dvy p"^ p{p Au Av) (2.3) 
Jo Jo Jo Jo 

where Au = Ux — Uy, Av = Vx — v. 



y-j i—^^ ^x ^y- 

a ra pUx rVx 



(Sx) = pV-2 dux dvx I dAu I dAv p^ p{p Au Av) 
'o Jo Jo Jo 

pV — 2 / dUx I dVx [integrand 1]^ 



^0 



JAd=0 



J Au=0 



where 



integrand 1 = — ^ (1 — pAuAv) exp{—pAuAv) 



Hence 



{Sx) = 1 - exp(-p a^) = l- exp{-pV) . (2.4) 

As p — 7- oo, {Sx) 1 and we write {Sx) ~ 1 to denote this. 

Consider now splitting up the interval X into four smaller intervals Ij, i = 1, ... 4, 



as shown in Fig. [Taj When computing the expected value of Sx one can split the 
integral up into the means of the actions of the four subintervals plus the "bilocal" 
contributions when x and y lie in two different subintervals. More concretely, given 
any subcausets, A and i? of a causal set C, we define the bilocal action, 

S[C- A, B] = N{A, B) - 2iVi(A, B) + AN^iA, B) - 2N^{A, B) (2.5) 

where N{A,B) is the number of elements m. Ar\ B and Nm{A,B) is the number of 
inclusive order intervals in C of cardinality m + 1 with top element in A and bottom 
element in B. Now let X and Y be submanifolds of spacetime M.. We define the ran- 
dom variable, ^m\X,Yi the Discrete Bilocal Action, via the sprinkling process: sprinkle 
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into M. at density p to obtain causet C with subcauset A{B) being that sprinkled 
into X(Y). For that reahsation, Sm;X,y takes the value S[C; A, B]. 
Note that Sm;X,x = Sx if X is a causally convex subset of A^j^ 
Now, consider X and its subintervals. If we adopt Sjj as simplified notation for the 
bifocal action Sx-x i , then we have 



{Sx) = J2{Sx^) + J2(%)- 



(2.6) 



i=l 



«J=1 
j<i 



The bilocal summands can be computed using the integral in Eq. (2.3) and adjusting 




(a) a and c {b and d) are the v- 
coordinate (it-coordinate) lengths 
of the sides of the subintervals 




(b) 



Figure 1: Splitting up a causal interval in 2D Minkowski to compute the action 



the boundaries. This yields 
(S21) = -2 



pb+d pb 

dvx / dUx / dvy / duy p"^ p{p Au Av) 

Jb Jo Jo 



-2 / dv, 

'0 Jb 



du^ 



[integrand 1] 



Au=Ux—b 



Av=0 



-1 + exp(— a6p) + exp(— adp) — exp(— a (6 + d) p) 



-1 



(2.7) 



^To simplify notation, we don't make the dependence on the density explicit. 

causally convex region, X, of A4 is one such that x,y Q X implies that the causal interval in 
Ai between x and y is a subset of X. 
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and 



(841) 



ra+c 
-2 / dVr 



I pa pb 

dUx / dvy / duy p'^ p{p Au Av) 
Jo Jo 



pa+c pb+d 

—2 / dvx I du^ 

J a J b 

1 — exp(— (a + c) (6 + d) p) 



integrand 1]^^^ 



- Av=Vx—a 



+ exp(— a (6 + d) p) + exp(— c {b + li) p)) 

+ exp(— (a + c)b p) + exp(— (a + c) dp) 

— exp(— a6p) — exp(— adp) — exp(— c6p) — exp(— cdp) 



(2.8) 



The three other bilocal contributions (Sjj) can be obtained from (S21) by changing the 



parameters appropriately. Putting together all parts of Eq. (2.6) one exactly recovers 



Eq. (2.1) 



Now, one can continue this game and split up the interval even further as in Fig. 
lb To compute the mean of the action one must again calculate 

9 9 

(Sx) = X](Sx.) + X](S.,). (2.9) 



i=l 



«J=1 
j<i 



We already know the contributions (S/-) ~ 1 and the bilocal contributions from two 
intervals that either share an edge or lie above and below a shared vertex {e.g. (S21) and 



S51) in Fig. lb ). It remains to compute the bilocal contributions from pairs of intervals 



such as (4,1), (7,1) and (9,1) in Fig. lb It turns out they consist only of exponential 



terms that are small when intervening intervals are large on the discreteness scale. In 
the limit of large density, we are left with a contribution of 1 for every subinterval, —1 
for every edge and 1 for every vertex. One could write 



{S)^F-E + V 



(2.10) 



where F denotes the number of faces i. e. intervals, E the number of edges and V the 
number of vertices. F—E+V is the formula for the Euler character of a polyhedron and 
motivates the question: Is the expected action (to some extent) a topological invariant? 
It is obvious that the formula can be applied to arbitrary causally convex regions of 
that can be tiled by causal intervals as long as each interval is large enough for the 
corrections to be negligible. It is not hard to verify that any such region will have a 
mean Discrete Action (S) ~ 1. So for example the region shown in Fig. 2a will give 
(S) ~ 1 but the region in Fig. 2b will not. 
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(a) causally convex (b) not causally convex 

Figure 2: Different regions constructed from causal intervals in 



3 Causally convex regions in M' 



The boundary of a causally convex region of can be spacelike in parts, but never 
timelike. If the region's boundary comprises straight line segments, such as the hexagon 



shown in Fig. 3b , then it can be divided up by null lines into a collection of intervals 



and causally convex triangles such as Fig. 3a Then the formula (2.10) will apply if 



the mean of the Discrete Action for a causally convex triangle tends to 1 in the infinite 
density limit. 



(a) causally convex triangle (b) causally convex hexagon 

Figure 3: Causally convex regions with boundaries formed from null and spacelike line 
segments 

First note that by Poincare invariance we can choose coordinates so that the spacelike 
edge of the triangle is at t =constant, and the apex lies at the origin. 
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Using null coordinates, as before, we have 



pL pvx rvx rux 

(Sa) = pV — 2p^ / dvx / dux / dvy / duyp{pAuAv) (3.1) 
Jo Jo Jo Jo 

where L = \/2V and V is the area of the triangle. This gives 

(Sa) = 1 + ^ + 0((pV)-^)^1. (3.2) 

We see that the mean DA of the triangle does indeed tend to 1 as p — )■ oo, though the 
corrections are not exponentially small. 

Now, consider a general causally convex region with a boundary whose spacelike 
portion is curved. So long as the discreteness scale is small enough - small compared 
to the radius of curvature of the boundary - we can tile the region with intervals and 
with causally convex approximate triangles along the spacelike boundary, all of which 



are large enough compared to the discreteness scale for the Formula (2.10) to hold 
approximately. We conclude that the mean of the DA for any causally convex region 
of will tend to 1 in the limit of infinite density. 

Is causal convexity necessary for the mean of the DA to be approximately 1? When 
a region, R C M^, is not causally convex, there will exist pairs of points x,y & R 
such that the causal interval in R between x and y is smaller than the causal interval 
between them in (the "diamond"). Since it is the volume of the causal interval in 
R which appears in the expression for the mean of the DA, one might expect this to 
disrupt the result and indeed it does. 

Consider the Discrete Action, Sn of a rectangle with edges parallel to the t and 
X axes. Analytic computation of the expectation value (Sn) is hard exactly because 
of the lack of causal convexity: the integral (1.3) breaks up into several subintegrals 
depending on the positions of x and y relative to the boundary. Therefore we use 
simulations to estimate the value. A sprinkling into a rectangle has three independent 
parameters that fully characterise the problem. One choice is the spatial width w, 
the height along the time-axis h and the sprinkling density The expectation value 
(Sn,^,h,p) niust be invariant under rescaling 

w — > X ■ w 

h ^ X-h (3.3) 
p ^ A-2.p. 



Fig. |4] shows simulation data for two different setups with power-law fits. Fig. 4a 
''Width w, height h and expected number of sprinkled elements N would be another choice. 
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100 200 400 800 

P 



(a) Simulation data for the action of a rectangle in for 
w ~ h — 1, varying density p with a power-law fit. Data 
averaged over 10^ to 10'' runs. Fit function: p°- ■ b. 




height h 



(b) Simulation data for the action of a rectangle in for 
w — 1, p — 100, varying height h with a power-law fit. Data 
averaged over 10^ to 10^ runs. Fit function: /i° • b. 

Figure 4: Numerical results for the action of a rectangle in M^. 



shows (S) for constant w and h and varying p, Fig. 4b for constant w and p and for 
varying h. Given the small relative error bars the power-law fits look quite convincing 
and we will assume that (Sn) can, at least in the regime covered by the simulations, 
be written in the form 

(Sn) = const ■ h'^w^p"'. (3.4) 



The scale invariance (3.3) demands a + j3 — 2'y = 0. From simulation 1 (Fig. 4a) one 
is tempted to deduce 7 = 1/2 and from simulation 2 (Fig. 4b) that a = 1. It follows 
/3 = 0. 

The fact that for constant p the width does not affect the value of the action whereas 
(Sn) oc h suggests that in general (Sn) contains boundary terms from timelike bound- 
aries only. We return to this question in the discussion section. 
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4 The flat cylinder 



In order to apply formula Eq. (2.10) to a causal interval, T^, of height T on a cylin- 
der with circumference L with L < T < 2L one might come up with a tiling into 
subintervals, Xj, i = 1, . . . 8, as shown in Fig. [5j Taking into account the topological 




Figure 5: Tiling of the interval Xc with L <T < 2L. a and b are the u and v coordinate 
lengths of the sides of the subintervals shown. 



identification, we have F = 8, E = 12,V = A thus yielding a predicted high-density 
expectation value of (Sx,.) ~ 0. However we have not shown yet that formula Eq. 
(2.10) is applicable to the cylinder. The division of the causal interval in Fig. Islhas 



been chosen such that formula Eq. (2.1) for the faces and formulae Eq. (2.7) and ( |2.8[ ) 
for the bilocal contributions of two intervals that share an edge or lie above and below 
a vertex can still be applied as the cylinder topology does not affect these cases. But 
the computation of contributions like (5,1), (6,1) and (8,1) differs from the Minkowski 
setup due to the nontrivial topology. 
Recall 

(Sx.) =pV-2K (4.1) 

where 

K = p^ [ (fyf d\ pipV^y) . (4.2) 

JXc JXcDJ+iy) 

In general, K can be split into a sum of terms, K = Ka depending on how many 
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homotopy classes of causal curves there are from y to x: 



Ka := p' / d'y / d'x p{pV,y) (4.3) 



where 



Jain) '■= & J^{y) \^ exactly a homotopy classes of causal curves from y to x} . 

(4.4) 

This split is motivated by the fact that V^y strongly depends on the number of homotopy 
classes of causal paths between x and y. For our interval, Ka = for a > 3. 
From Fig. [5] we see the relation between a, b, T and L is: 



a = {T-L)/V2 

b = {2L-T)/V2 (4.5) 



The causal volume V^y for x G Jaiu) for a > 3 is at least (a + 6)^ so is suppressed 
by at least exp(— p(a + 6)^) and can thus be neglected as L = \/2{a + b) is assumed to 
be large in discreteness units of p~5. 
The values for Ki and K2 are [12] 

Ki = + - exp(— pa^) — (1 + pab) exp(— pa(a + b)) + corr. 
2 

K2 = -7 z-:^ +exp(-pa(a + b))[l + pab 

[a + Oj^p 

+ 7 ^-r^ ((6 + 2p(a + 6)(2a + b) - p\a + bYb"" + p='(a + bfab"") 

(a + bpp'^ ^ 

— 2exp(— pa(a + &))(3 + 4pa(a + 6) + 2p^a^{a + 6)^))] + corr. (4.6) 

where "+ corr." stands for neglected terms suppressed by exp(—p(a + 6)^). However we 
will keep terms with factors exp(— pa^) and exp(— pa (a + b)) since for T only slightly 
larger than L the value of a will be very small and these terms are then significant. 
The overall action is 

(Sj^) = — exp(— pa^) + 2(1 + pab) exp(— pa(a + b)) 
4 

+ 7 7V7^ + exp(-pa(a + b)) [I + pab 

[a + byp 

+ 7 \rr^ ((6 + 2p(a + b){2a + b) - p\a + bfb^ + p'^ia + bfab'') 

[a + oj*p^ 

-2exp(-pa(a + 6))(3 + 4pa(a + 6) + 2p^a^(a + 6)^))] + corr. . (4.7) 
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For T > 2L consider a division of the interval into regions 1 and 2 as shown in Fig. 
[6) The expected DA is the sum of the expected actions for regions 1 and 2 and the 
bilocal contribution (S21). 



r 




2 



> 2L 





J 



1 



V 




Figure 6: Division of interval when T > 2L. 



It can be shown fri\ that the expected action for region 1 and the bilocal contribution 
cancel (up to exponentially small terms) and the result is just given by the expected 



to simulation results. For T ^ L the action approaches the Minkowskian limit 1. For 
T only slightly greater than L the exponential terms dominate and cause a downwards 
spike. As T — 7- 2L the non-exponential correction, ^ (which comes from K2) dom- 
inates. However this also tends to zero in the limit p — )■ 00 so (Sx^) ~ as initially 
predicted. Indeed it can be shown explicitly that the bilocal contributions from pairs 
of intervals that do not share and edge or vertex tend to zero as p — ?■ 00 and so the 
formula F — E + V can be applied to intervals of the cylinder. More generally, the 
previous argument regarding null tilings of causally convex regions of can be given 
here, and we conclude that (S) ~ for general topologically non-trivial causally convex 
regions of the cylinder. 



action of region 2 which can be obtained from Eq.(4.7) by setting a = L/\/2,b = 
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-0.5 



T 



T 



simulation, N=100 
simulation, N=200 
analytic, N=10{) 
analytic, N=20{) 




Figure 7: The expected action of a cylinder-interval for L 
200 compared with simulation results. 



1, (iV) = 100 and (A^) 



5 The flat trousers 



We investigate now a causally convex neighbourhood of the fiat 1+1 trousers spacetime 
in which two S'^'s join to form a single S^. The trousers spacetime is a piece of with 
cuts and identifications as shown in Fig. |8| Although the singularity, P, at which the 
topology changes is by some definitions not strictly in the spacetime since the metric 
degenerates there, nevertheless the causal order is well defined at the singularity: it is 
clear what the causal past and causal future of P are. Therefore we will consider P as 
a point of the manifold. Note that in any sprinkling into the trousers almost surely no 
element will be sprinkled at P. 



► ♦ 4' 



Figure 8: The trousers spacetime. P is the singularity - all three instances of P are 
identified - and the shaded region is a neighbourhood of P. There is a vertical cut 
down from the central copy of P with the two legs identified as shown. 



Let M denote the neighbourhood of P shown as the shaded region in Fig. [8j It 
consists of two fiat intervals each with P as their midpoint, identified across "branch 
cuts" from P to their past tips. N is topologically a disc if P is in included the manifold 
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and if the formula (2.10) holds then the expected DA of M would be equal to 1 in the 
limit of large density. 

Let the volume (area) of each of the two intervals be Aa? and consider the null tiling 
into 8 intervals, Xj, i = 1, . . . , 8, shown in Fig. [9] The interval Xi comprises the two 
triangles labelled 1' and 1" and the interval X2 comprises the triangles labelled 2' and 
2". Adopting the same notation for the bilocal discrete action of two intervals used in 



(2.6) we have 



(SAr) = I](Sx.) + ^(S.,). (5.1) 



i=l *ii=l 




Figure 9: Null tiling of Af into 8 intervals. 

For each i, (Sx^) ~ 1. The bilocal terms are nonzero when the intervals Xj and Xj 
share an edge and in that case (Sy) ~ —1. There are 8 edges so these contributions 
cancel the contributions of the 8 individual intervals. The only other nonzero bilocal 
terms are (Sjj) where i = 5,8 and j = 1,2 and their sum is the contribution of the vertex 
at the singularity. These 4 terms are equal by symmetry so we have (S_^) = 4(851). 



The causal interval between x G X5 and y G Xi is shown in Fig. [TO] and we deduce 
that 



/'2a p2a pa pa 

(S51) = -2 / dux I dv^ / duy / dvy p'^ p{pVxy) . (5.2) 

Ja Ja Jo Jo 

Vxy = AuAv — {v^ — a) (a — Uy) . (5.3) 



where 



This gives 



(Sat) = 4 ln{pa') + 4(7 - 1) + O ( ) (5.4) 
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Figure 10: The causal interval between a; G X5 and ?/ G Xi is depicted in grey. 



where 7 is Euler's constant. We see that the expected DA of the neighbourhood of 
the singularity does not tend to 1 or any constant but grows logarithmically with the 
density. 

6 Discussion 

We have shown that in the limit of infinite density, the mean of the Discrete Action 
will be 1 for any causally convex region of including regions whose past and/or 
future boundaries contain spacelike segments. Since these spacelike segments may 
have nonzero geodesic curvature, the constancy of the mean of the DA suggests that 
it contains no contribution from the past or future boundaries. 

Indeed a handwaving argument can be given as to why this should be so, even when 
Ai is curved. The boundary of a causally convex region U <Z M. consists of a future 
boundary and a past boundary intersecting in a co-dimension 2 spacelike surface. There 
is no timelike portion of the boundary. The mean of the DA is a double integral over 
U which can be done in either order. The integrand is a retarded 2-point function. 



where ?/) = 1 if ?/ G J^(x) and otherwise. Let us assume the density is high 
enough that a sprinkled causal set can capture the curvature of A^, i.e. at each point 
y E Ai there is a local inertial frame in which the curvature components are small 
compared to the density. If we do the x integration first, at fixed ?/, then the resulting 
function pL{y) is approximately \R{y) unless y is too close to the future boundary. If 
it is within length p^^ of the boundary then the range of the x integration will not 



pC{x,y) 



P 



6'~'\x,y)-2p'p{pV,y)C{x,y) 



(6.1) 
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be large enough for the approximation to hold [1]. Then we do the integration over 
y to get approximately the usual Einstein-Hilbert bulk term together possibly with 
some contribution from the integral over the points y close to the future boundary, i.e. 
possibly some kind of future boundary term. But there is no contribution from the 
past boundary at all. Now reverse the order of integration: do y first and then x. Now 
there appears to be no contribution from the future boundary. This can only happen 
if neither boundary contributes. So the only points where some boundary contribution 
can come in, is from the points which are close to both past and future boundaries 
i. e. from the spacelike co-dimension 2 "corners" where the past and future boundaries 
intersect. The argument holds when the past and future boundaries are partly spacelike 
as well as when they are wholly null. There is no reason, from this argument, that 
timelike boundaries could not contribute however and we saw evidence that they do 
from the results for the rectangle. 

This heuristic reasoning would have to be backed up with further evidence from 
simulations of the Discrete Action but it suggests that there is no Gauss-Bonnet formula 
for the 2D Discrete Action. The 2D Gauss-Bonnet Theorem can hold because, as the 
geometry of the bulk surface is varied, the extrinsic curvature of the boundary changes 
and the right combination of bulk and boundary terms can remain constant. In 2D the 
co-dimension 2 "corner" is an , i.e. 2 points, and if the only boundary contributions 
are from these 2 points, these couldn't compensate for the changing bulk term. Another 
reason not to expect the DA to satisfy a Gauss-Bonnet formula is that it appears that 
the appropriate Lorentzian analogue of the Euclidean formula is of the form "bulk 
term + boundary term + corner terms" = ^Tiix rather than 2ttx [l3l [H] (see also 
[ini [16]). Both the bulk and boundary terms are real but the formula can hold because 
the corner contributions are Lorentzian angles which can be complex. However, the 
Discrete Action is real. 

This putative lack of boundary terms could explain why the expected DA for any 
causally convex region of is the same. The continuum bulk term is zero. If the mean 
of the DA is indeed close to the continuum bulk term plus only a contribution from the 
corners then that should be the same for all causally convex regions. Presumably, 
the difference for the neighbourhood of the singularity of the trousers comes from a 
boundary effect of the non-standard causal structure around the singularity which has 
a double lobed past and future. These issues all remain to be investigated. 

There are a large number of open questions. What does happen in 2D curved space- 
times? Will the results bear out the conjecture that the expected DA is approximately 
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the Einstein-Hilbert term plus a constant from the corner? What happens in higher 
dimensions? There are analogues of the 2D Discrete Action in 4D [1] and 3,5,6,7D and 
higher [5]. One would expect, for example, that the mean of the DA of an interval in 
M*^ would be proportional to the volume of the S"^^"^ corner. 

Although one need not take any position on quantum gravity to find interest in the 
Discrete Action as a random variable defined for a continuum spacetime - one need not 
consider the discreteness of the causal sets that arise in the definition of the Discrete 
Action to have any physical basis - its main application is likely to be in the causal set 
approach to quantum gravity. So, what is the significance in quantum gravity of the 
results for the interval, trousers and rectangle? For example, the result for the rectangle 
suggests that the expected DA contains boundary contributions proportional to the 
length of any timelike boundary. Can we use the DA to give an argument against the 
appearance of "holes" and "edges" in spacetime? Or for or against topology changing 
processes such as the trousers? 

A major open question is how the fluctuations in the DA behave as the density gets 
large: we should stress that the results reported here are all concerning the mean of 
the DA. For a typical sprinkled causet, how far is the DA from the mean? Preliminary 
results show the fluctuations grow as the density gets large [12], contrary to the hope 
expressed in |lj and this needs to be studied further. To tame the fiuctations it may be 
necessary to introduce a mesoscale between the discreteness scale and the observation 
scale [171 H]. Further work is needed to illuminate these issues. 
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